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The electrophoretic motion of a conducting particle, driven by an induced-charge mechanism, 
is analyzed. The dependence of the motion upon particle shape is embodied in four tensorial 
coefficients that relate the particle velocities to the externally-applied electric field. Several families 
of particle shapes, whose members are unaffected by the field, are identified via use of symmetry 
arguments. Other particles translate and/or rotate in response to the imposed field, even if their net 
electric charge vanishes. The coefficients are represented as surface integrals of the electric potential 
over the particle boundary, thereby eliminating the need to solve the flow field. 



Traditionally, the term electrophoresis applies to the 
motion of colloidal particles through an electrolyte solu- 
tion. This motion results from the interaction between an 
applied electric field, say Eoo , and the net ionic charge ac- 
commodated in the Debye cloud surrounding the particle 
boundary (this charge being equal in magnitude and op- 
posite in sign to that adsorbed at the particle surface). In 
many situations, the Debye-layer thickness is small com- 
pared with particle size. The electrokinetic flow is then 
described by an approximated equation set, obtained via 
use of singular perturbations []J. In that asymptotic de- 
scription, the Debye screening length assumes zero size. 
Thus, the electric potential ip and velocity field v in 
the electrically-neutral fluid domain are respectively gov- 
erned by Laplace's and Stokes' equations, whereas the 
presence of of the layer is reflected by the boundary con- 
ditions that apply at the particle boundary S. These con- 
sist of the no- flux condition, n (x s ) -E (x s ) = 0, governing 
the normal component of the electric field, and the slip 
condition, vs (x s ) = — (e/n) £ (x s ) E (x s ), governing the 
relative fluid-particle motion. These conditions respec- 
tively reflect the particle impermeability to ionic current 
and the net action of the electric field on the charged 
Debye layer. Here, n is an outward-pointing unit vector 
normal to S, x s is a position vector located on S, e and 
fi respectively denote the permittivity and viscosity of 
the solution, the "zeta potential" ( constitutes the local 
value of the potential jump across the Debye layer, and 
E = — V<p is the electric field. 

The preceding description is adequate for dielectric me- 
dia, capable of carrying bound electric charge. Electroki- 
netic flow can also be generated about a conducting body; 
such flow, however, is driven by a polarization mecha- 
nism, made possible by the charge mobility. As in the 
dielectric case, the electric charge is also situated on the 
surface; this, however, is not the result of any physic- 
ochemical bonding, but rather a simple consequence of 
Gauss law 2]. Moreover, the surface charge in a conduc- 
tor distributes over its surface so as to guarantee zero 
electric field in its interior. The resulting dependence 
of the zeta potential upon E^, gives rise to a nonlin- 
ear slip mechanism, which is responsible for a rich va- 
riety of physical phenomena. The most notable feature 



in this induced electrokinetic effect is the generation of 
electrokinetic flow about an initially-uncharged particle, 
this flow having no counterpart in conventional linear 
electrophoresis. 

Flows driven by this induced mechanism, addressing 
both DC and AC external fields, have been discussed 
mainly in the Ukrainian literature (see 0). A detailed 
description of induced-charge electrokinetics appeared in 
a recent paper by Squires and Bazant Q. The au- 
thors, discussing a broader context of flows about polar- 
izable media, coined the terms "induced-charge electro- 
osmosis" (ICEO) and "induced-charge electrophoresis" 
(ICEP) to respectively describe the associated fluid and 
particle motion. They also suggested |j| that the com- 
bination of nonlinear slip with symmetry breaking can 
lead to useful configurations which can perform a variety 
of microfluidic operations (such as pumping and mixing) . 
Their suggestions focus upon electro-osmotic flows rela- 
tive to fixed conducting elements (which may be held at 
a constant electric potential) , rather then electrophoretic 
motion of freely-suspended particles, which is addressed 
herein. Nevertheless, electrophoretic analyses tend to 
provides insight to the inverse problem of electro-osmosis. 

The simplest ICEO situation is that of an initially un- 
charged conducting particle which is positioned in a uni- 
form and constant applied field E^ . During a short time 
interval, a Faradaic current charges the Debye layer in 
proximity to the particle surface, with an equal and op- 
posite charge being set up on the surface itself. The 
induced potential difference across the layer, namely the 
zeta potential, is proportional to E^. Accordingly, the 
electrokinetic slip that occurs at S is quadratic in E^. 
Given the linearity of the flow problem, so are all the 
variables describing the resulting fluid (and perhaps par- 
ticle) motion. 

The prototypical ICEO problem, which was investi- 
gated in 0| as a model for the electrokinetic flow about 
mercury drops, entails a spherical particle. Owing to the 
fore-aft symmetry of this flow 6] , no hydrodynamic force 
(or torque) is exerted on the particle. This dipolar sym- 
metry is obviously absent in flows about more general 
particle shapes. In principle, any deviation from a spher- 
ical shape may lead to hydrodynamic force and torque, 
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with a consequent particle motion. Indeed, such motion 
has been experimentally observed Q. 

The purpose of this paper is to elucidate the inter- 
esting properties of this nonlinear phenomenon. I con- 
sider here a relatively simple case, addressing the elec- 
trophoretic motion of a freely-suspended conducting par- 
ticle in a uniform and constant externally-applied field. 
(The present results are easily generalized to the case 
of AC fields.) The analysis begins with the case of an 
initially-uncharged particle and is followed by discussion 
of non-zero net charge effects. 

Consider an uncharged conducting particle of an ar- 
bitrary shape and of linear dimension a, which is posi- 
tioned in an electrolyte solution in the presence of an 
externally- applied constant electric field, Eqo — K30E (E 
being a unit vector in the direction of Eqo). When de- 
scribing the governing field equations, the position vector 
is conveniently measured from an (arbitrary) point O at- 
tached to the particle. It is also convenient to describe 
the particle motion by the instantaneous velocity of point 
O, Uo, together with its angular velocity f2. 

For conducting particles which are impermeable to the 
current carrying ions (as is usually the case) the no-flux 
condition, n ■ V<^ = 0, still holds. Thus, the electro- 
static problem governing ip consists of this condition, 
together with Laplace's equation in the fluid domain, 
V 2 ip = 0, and the far field condition, V</5 — » — Eqo- This 
problem (which also describes the field about a compa- 
rable dielectric particle) is linear and homogeneous in ip 
and Eqo, from which it follows that ip must be a lin- 
ear function of Eqo- The attendant flow field is gov- 
erned by the Stokes equations, the slip condition at <S, 
v = Uo + ft x x + (s(/fi) V</3, and the attenuation con- 
dition at large distances from S, v — > 0. For a freely 
suspended particle, the requirement of vanishing hydro- 
dynamic force and torque completely determines the flow 
field, as well as the values of Uo and fl. 

As in common in the literature (see Q), it is assumed 
that the zeta potential is sufficiently low to be propor- 
tional to the local surface charge density q, 



Here, A is the Debye-layer thickness of the Debye layer. 
With no loss of generality, the uniform value of the elec- 
tric potential in the particle interior is chosen as zero. 
Accordingly, the zeta potential is related to the potential 
value just outside the layer, ( = — ip\g. It is the nonlin- 
ear slip term (pVtp which is responsible for the peculiar 
characteristics of ICEP. Since the slip is the driver of the 
flow, this term results in a quadratic dependence of v 
upon Eqo. Owing to the linearity of the flow problem, 
the particle translational and rotational velocities must 




FIG. 1: An axisymmetric particle possessing fore-aft symme- 
try: in the absence of net charge, the particle remains station- 
ary; in the presence of an initial charge, the particle translates 
but does not rotate. 

possess the following structure: 

Ed 

Uo = — Co : EooEoo, (2a) 
A 4 

n = — D : E^E^. (2b) 

Here, Co is a dimensionless third-order tensor (which 
depends upon the position of O) and D is a dimension- 
less third-order pseudo-tensor. These objects are intrin- 
sic geometric properties of the particle, independent of 
its size and orientation relative to the imposed field. In 
a sense, they constitute a lumped description of the elec- 
trophoretic motion. 

It is actually possible to demonstrate a variety of sym- 
metry properties, expressed through the coefficients Co 
and D, for various families of particle shapes. These 
properties may be obtained without the need to explic- 
itly solve the governing field equations. The simplest ex- 
ample involves isotropic particles (such as spheres). For 
that family, no vector (or tensor) can specify the parti- 
cle geometry. Since no isotropic third-order tensor ex- 
ists it is impossible to construct a candidate for Co- 
Moreover, the only isotropic third-order pseudo-tensor is 
e, the alternating triadic, which yields zero when con- 
tracted with the symmetric dyadic EooEoo. Accordingly, 
isotropic particles neither translate nor rotate. 

A larger family is that of all axisymmetric particles 
possessing fore-aft symmetries, such as spheroids and fi- 
nite cylinders of circular cross section (see Fig. 1). These 
particles are characterized by a single vector e attached 
to their axis of symmetry. The only third-order tensors 
which may be constructed using e are eee and various 
permutations of le (I being the idem-factor). These ten- 
sors, however, are odd functions of e, and do not satisfy 
the invariance to reversal of e required from a body pos- 
sessing a fore-aft symmetry. Moreover, it is impossible 
to construct a third-order pseudo-tensor using the single 
vector e. Accordingly, such particles behave like spheres: 
they neither translate nor rotate. 

It appears that the generation of particle motion re- 
quires a stronger symmetry breaking. Consider the case 
of axisymmetric particles lacking fore-aft symmetry, such 
as cones and hemispheres. Again, these particles are 
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characterized by a single vector e attached to their axis 
of symmetry, whence the only candidate for D is again 
e, which means that such particles do not rotate. Both 
XJo and Co are therefore independent of O, and may be 
respectively denoted by U and C. As with the previous 
family, the only candidates for C are eee and various per- 
mutations of le. For bodies lacking fore-aft symmetry, C 
is not necessarily invariant to reflection. Thus, such par- 
ticles may translate (without rotation) under the action 
of an external field. 

In principle, both Co and D are obtainable from the de- 
tailed knowledge of the velocity field. In what follows, in- 
tegral expressions for these tensors are derived in terms of 
the electric potential distribution. These expressions are 
calculated using the Lorentz reciprocal theorem , thus 
avoiding the formidable calculation of the flow problem. 
For simplicity of presentation, this procedure is demon- 
strated for non-rotating (e.g., axisymmetric) particles, 
for which the angular momentum balance is automati- 
cally satisfied (the generalization to arbitrary particles 
is trivial). In what follows, it is convenient to employ 
a dimensionless notation, wherein length variables are 
normalized with a, velocities with eaE^/fi, forces with 
ea?E 2 <l , stresses with eE^, and the electric potential with 
aEoa. 

Using this notation, the slip condition on S appears 
as v = U — (p"V(p (~ denotes here the dimensionless 
counterpart of the relevant physical quantity) . The lin- 
ear flow problem is decomposed into two separate parts. 
The first describes pure translation of the particle with 
a velocity U (in the absence of any slip); in the second 
the particle is held fixed, and the flow is generated by the 
imposed slip on S, v = — <f> V£>. In both problems the ve- 
locity field satisfies the Stokes equations and attenuates 
at large distances from the particle. 

The force resulting from the first flow field is given by a 
Stokes drag expression, Fi = — K • U, where the positive- 
definite dimensionless translation tensor K is an intrinsic 
property of the particle geometry . The force resulting 
from the second flow field is provided by the following 
quadrature: 

F„ = - j> dAh- S f • Vlp (3) 

which follows from a variant of the reciprocal theorem 
[To). Here, dA is a dimensionless area element (nor- 
malized with a 2 ) and the (dimensionless) "translational" 
stress triadic field £ (defined in yj) depends only upon 
the particle geometry and the position vector of the fluid 
point on S: explicitly, Syfc denotes the ij component of 
the stress tensor that would result from a pure transla- 
tion of the particle with a unit velocity in the k direction. 
(The superscript t signifies right transposition.) 

The requirement of a force-free particle, Fi + Fn = 0, 



yields the following result: 

U = -K _1 • <j> dAh ■ • ipVip. (4) 

Owing to the linearity of the electrostatic problem, it is 
clear that (p = <I? (x) • E, where 3>, a dimensionless vec- 
tor function of position, can only depends upon particle 
shape. Use of the definition (12a[) furnishes the requisite 
representation, 

C = -K _1 • <j> din • • V3> (5) 

It is emphasize that this result holds only for a particle 
which does not possess an initial charge. 

Consider now a particle possessing an initial charge Q 
(or, equivalently, a net charge Q). In this situation, the 
zeta potential is not quadratic in Eoo and the relations 
(0) are invalidated. Define an "effective" potential, ip = f 
QX/Ae, where A = <f> s dA is the particle surface area. 
While it is customary to identify i\> with the zeta potential 
[cf. l|T]l]. I here avoid any such a priori interpretation. 

The flow about an initially-charged particle under the 
action of an applied field has been discussed in 0. 
Given the linear-capacitor approximation additivity 
in charge implies additivity in the zeta potential. Accord- 
ingly, the initially-charged situation was represented as a 
combination of two different flows: the first is associated 
with a quadratic slip condition, and the second is that 
involved in a linear electrophoretic motion correspond- 
ing to a uniform zeta potential £ = ip. It is important to 
note, however, that the situation discussed in [3j — with 
a view towards microfluidic devices — is focused upon 
the specific case of a circular cylinder (a similar repre- 
sentation is also valid for the case of spheres). Is that 
interpretation valid for more general shapes? 

The additional scale ip modifies the translational and 
rotational relations (J2J , which, following dimensional ar- 
guments, are replaced by: 

SOj s s 

Uo = — A : E^Eoo + -VM • E^ + — ^ 2 Po(,6a) 
/i /i /ia 

n = -B : EooEoo + — V>N • Eoo + -Ar ^R- (6b) 
/i /ia /ia 

The dimensionless objects Aq, Mo, and Po respectively 
denote a triadic, dyadic, and a vector, whereas the di- 
mensionless objects B, N, and R respectively denote a 
pseudo-triadic, pseudo-dyadic, and a pseudo- vector. All 
six quantities are intrinsic geometric properties of the 
particle, independent of the applied field and the value 
of ip. 

In the case Eoo = the translational and rotational ve- 
locities are respectively proportional to Po and R. Since 
the applied field is the driver of the motion, Po and R 
must vanish in that case. These quantities, however, are 
independent of Eoo, and must therefore vanish identi- 
cally. Similar logic can be applied via the case ip = 0, 
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which yields relations similar to J5J. Since Ao and B are 
independent of ij), they must be respectively identical to 
the coefficients Co and D, which describe the motion of 
particles possessing zero net charge. 

The role of Mo and N is elucidated via the asymptotic 
limit aEoz/^p <g; 1, which corresponds to the quadratic 
effect (associated with the action of Eoo on charge in- 
duced by Eqo) being negligible relative to the linear one 
(associated with the action of Eqo on the initial charge 
distribution). In this case H6a(l - (|6b|l degenerate to: 

U = -V^Mo-Eoo, (7a) 
A* 

fl = — V>N • Eoo- (7b) 

These are identical to the mobility relations governing 
conventional electrophoresis (with no induced charge ef- 
fects), except for an important difference which will be 
identified shortly. Since Mo and N are independent of i/j 
and Eoo (and , specifically, of the above-mentioned limit), 
they may be obtained using a linear electrophoretic anal- 
ysis. It is therefore tempting to say (see pf) that the 
mobility relations J7J) are the same as those for a lin- 
ear electrophoresis of particles having a zeta potential 
i/>, which would imply that Mo and N are governed by 
the Smoluchowski relation U = (eip / ijl)'Ei 00 and $1 = 0. 
That would yield M = I and N = 0. 

It is important, however, to emphasize that the linear 
electrophoretic motion in the limit aE^/ip <C 1 corre- 
sponds to a charge distribution in a conducting particle 
in the absence of any imposed field. Only for the case of 
a sphere (or a circular cylinder in two-dimensional flows) 
this distribution is uniform, corresponding to a uniform 
value (namely ijj) of the zeta potential. Thus, Mo and N 
would not be identical in general to those predicted by 
the Smoluchowski relations. (Indeed, linear electrophore- 
sis of non-uniformly charged dielectric particles is 
qualitatively different than that observed in the uniform 
case.) It appears that the traditional focus on spherical 
and cylindrical bodies leads to a loss of physical effects 
which characterize nonlinear slip. As such, conclusions 
drawn from the analyses of such specific shapes are rather 
restrictive and non-representative. 

As in the case of particles possessing no net charge, 
for which the relations J2Jl apply, it is possible to deduce 
symmetry properties for the motion (namely the values 
of Mo and N) of initially-charged particles. With ax- 
isymmctric particles, for example, Mo must be a linear 
combination of I and ee, whereas N must be proportional 
to e-e. In the low-field limit, aE^/ip <C 1, such particles 
do not rotate, but in general may translate. 

It is a simple matter to represent both Mo and N as 
quadratures of the electric field, without the need to solve 
the flow problem. As with the preceding analysis, which 
has led to the representation JSJ), the present derivation 
focuses upon non-rotating (e.g., axisymmetric) particles, 



for which N = and Mo is independent of O (again, 
the generalization to the general case is trivial). In the 
limit aE^/ip <C 1 the slip velocity is given by (s(/fi) V<p, 
with £ being the zeta potential associated with the ini- 
tial charge distribution (which is a function of x s ). It 
is convenient to use the previously employed dimension- 
less notation, with which thc_ present slip_ condition ap- 
pears as v = U + (ip/aE^CVft (with C = C/V0- The 
concomitant flow field is decomposed into two parts, the 
first describing pure translation with a velocity U, and 
the second describing a stationary particle on which the 
velocity slips with the prescribed value (ip/aE^ ) (Vtp. 
A similar procedure to that leading to 10} then yields 

U= -l-K" 1 - <fdA(h-^- V£. (8) 

Use of the linear representation, <p = <& • E, together with 
the definition (|7a.|> furnishes the desired expression, 

M = K _1 • <j> dA(h ■ • V*. (9) 

It is re-emphasized that the zeta-potential distribution 
appearing in @ corresponds to a surface charge arrange- 
ment on a conducting particle in a vacuum (in the ab- 
sence of any applied field), whereas the electric field is 
identical to that about a dielectric particle in a vacuum. 

For a sphere of radius a, the initial charge spreads uni- 
formly, giving rise to a charge density Q/A. Thus, ip is 
identical to the zeta potential, and £ = 1. The potential 
distribution under the action of an applied external field 
Eqo is given by a combination of a uniform field and a 
dipole, (p = — (l + a 3 /2r 3 ) x • Eoo, from which it follows 
that V*| 5 = (3/2) (fin - I). Also, from the solution 
of the Stokes equations for a translating sphere 0, it is 
known that n • S = (3/2) I and K = 6ir\. Use of the 
integral relations § s dA — 4ir and / 5 cL4nn = (Att/3) I 
readily yields the Smoluchowski relation M = I. 
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